Lognormal properties of SGR 1806-20 
and the Possibility of a Quiescent 
Population of other SGR sources 
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CN ! Abstract. 

Monte Carlo simulations of long SGR event sequences based on lognormal distribu- 
■ tions with a range of time intervals and intensity distribution parameters have been 

investigated. The main conclusions are that the majority of SGRs with properties 
similar to SGR 1806-20 have been detected but SGRs with mean waiting times much 
longer than SGR 1806-20 remain to be discovered. A large decrease in the probability 
for detection of an SGR source results from a relatively small increase in the distri- 
5o ! bution parameters obtained for SGR 1806-20. A new breed of experiments with very 

&\ • long observation times are required to search for this type of source. 
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The lognormal properties of the soft repeater SGR1806-20 have been previously 
reported [2,3,1]. In particular, both the time interval between repeater events and 
the luminosity function of the source were fit with lognormal distributions [4]. This 
analysis used the data-base of 111 events detected by the International Cometary 
Explorer (ICE) mission [5]. 

While the present number of events observed from the other three sources [6-8] 
does not allow any detailed analysis, the intervals between successive events of SGR 
0526-66 [9] is also suggestive of lognormal behaviour. 

The relationship between the number of active (i.e observable) sources and the 
true number of SGRs in the galaxy is one which is the subject of some debate 
[10,11]. If the time intervals between SGR events proves to be lognormal with a 
wide range of means and variances then there may be long quiescent periods where 
the source could be undetectable, leading to an underestimate of the population. 



Simulated SGR source activity over 500 year period. 
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FIGURE 1. Four separate Monte Carlo simulations of duration 500 years of SGR1806-20 activ- 
ity, generated from lognormal distributions for recurrence intervals and intensities with parameters 
as determined from SGR1 806-20 [1]. The detector model is described in the text. The long gaps 
in activity of the source are the contributions from the tail of this highly skewed distribution. 

SIMULATIONS 

Monte Carlo simulations of SGR event sequences with a variety of distribution 
parameters were generated. The simulations were performed using Octave 2.0.5 
for UNIX. Randomly generated standard normal variates, X, were transformed 
to lognormal variates using the relationship Y = e aX+ ^ where Y is lognormally 
distributed with parameters \i and a. The effect of the detector on the observations 
was simulated by rejecting events below a preset intensity level to mimic a threshold 
and by rejecting at random some predefined fraction of events to mimic a less than 
100% live time. 

The parameters of the two lognormal density function which were fit to the 
distribution of recurrence intervals and intensity for SGR1806-20 were respectively 
/iw = 9.64, aw = 3.44 and /ii = 3.3, o\ = 1.3 [1]. To investigate how the probability 
of observation with a given time interval depended on the the parameters for the 
waiting times (/iw an d cr w ) and the parameters for the event intensity (/xi and oi), 
two large simulations were performed. In the first, the waiting time parameters were 
varied over a small range relative to the SGR1806-20 values (while the intensity 
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FIGURE 2. This surface plot shows how the probability of observing an SGR source varies 
with the parameters of the waiting time distribution(^ w and aw)- The probability given is per 
0.5 years. The contour lines show the drop in probability with /i\y and o\y. The indicated cell 
contains the value for SGR 1806-20 with /iw = 9.6 and <tw = 3.4. 



distribution was fixed), forming a grid of 14 x 14 cells, each with a different value for 
/iw and <T\y but the same value for fi\ and o\. For each cell on the grid ten 500 year 
simulations were performed (as for the sequences produced in Figure 1) and the 
mean number of events per 0.5 year interval was determined. The second simulation 
was performed with varying fii and <7i(the intensity distribution parameters) but 
it followed the same procedure - the region of parameter space to be explored was 
divided into a grid and ten 500 year simulations were performed at each grid point 
(each grid point having a different /ii and o\ but the same /iw and ow)- The 
resulting distributions of probability of observation for each of the two grids are 
presented in Figures 2 and Figure 3. 



DISCUSSION 

Figure 2 presents a surface plot of the probability of detecting an SGR source 
within 0.5 years, given that the event intervals are distributed with a lognormal 
probability and the intensity distribution is fixed to that observed for SGR1806-20. 
Note that the parameters for SGR1806-20 (/zvk and aw) an estimated probability 
of 30% for observation in 0.5 years and that the probability falls quickly with 



FIGURE 3. This surface plot shows how the probability of observing an SGR source varies with 
the parameters of the intensity distribution(/iiand a{). The probability given is per 0.5 years. The 
indicated cell contains the value for SGR 1806-20 with /ii = 3.3 and <j\ = 1.3. 

increasing /^w and ow- For /z\y = 11 ~ 1.2/xJy 06 , o"w = 4.5 ~ 1.5cx^ 06 , which is 
equivalent to a mean of 1.5 x 10 9 s, the probability of observing the source in 0.5 
years has fallen to less than 5%. 

In contrast, Figure 3 shows a much flatter more uniform behaviour over a similar 
variation in the intensity distribution parameters. Note that the parameters for 
SGR1806-20 again give us an estimated probability of 30% for observation in the 
0.5 years of observation. The grainy nature of the surface is due to the much 
lower variation in probability across the range of intensities chosen (from 50% of 
SGR1806-20 to 300% for each distribution parameter). 

The lognormal distribution arises in statistical processes whose completion de- 
pend on a product of probabilities, arising from a combination of independent 
events [12]. Lognormal statistics have previously been used in connection with 
gamma-ray bursts [13,14]. In the context of this investigation the physical sig- 
nificance of this statistical behaviour may lie in the connection between SGRs 
and neutron stars where a similar statistical analysis was presented for the micro- 
glitches from the Vela pulsar. The time separation and the intensity of these small 
(| Av/v |~ 10 -9 ) frequency adjustments were both compatible with lognormal dis- 
tributions, and there was no correlation between waiting time and intensity, just 
as observed with SGR1806-20 [5]. This result, combined with the identification of 
X-ray point sources [15,16] embedded in plerion-powered SNR [17] as counterparts 



to the SGR sources, suggests structural adjustments in neutron stars may be the 
cause of SGRs. 

CONCLUSION 

The activity of sources with mean recurrence times similar to or much longer 
than SGR 1806-20 was investigated using Monte Carlo simulations. The results 
of the simulations indicate that there could exist a significant population of SGRs 
with mean waiting times considerably longer than SGR1806-20 that remain undis- 
covered. Structural adjustments in neutron stars may be responsible for this be- 
haviour. Finally, a new breed of experiments with very long observation times will 
be required to search for this type of source. 
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